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Abstract: The application of gravity anomalies for gravimetric geoid model determination has necessitated the computation of 
normal gravity on the WGS 84 ellipsoid using the International Gravity Formula. For local gravimetric geoid model determination, the 
gravity anomalies are computed on the local ellipsoid adopted for geodetic computation in the area/region of study. To determine 
precise local gravimetric geoid model in countries where Clarke 1880 ellipsoid is adopted for geodetic computation, the theoretical 
gravity must be computed on the adopted geodetic computation surface. As a result, this paper derives a local theoretical gravity model 
on the Clarke 1880 ellipsoid for practical local geoid model determination. The model was derived using the Clarke 1880 ellipsoid 
parameters, angular velocity, gravitational constant and the Somiglinana's closed formula for normal gravity. The derived model was 
presented in two forms, model A and model B. The two forms of the derived model were compared to determine their suitability as well 
as agreement. The comparison results show that the two forms of the model agree. Thus, any of the model forms can be applied for 
accurate normal gravity computation where Clarke 1880 ellipsoid is adopted as a computation surface. 
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INTRODUCTION 

One of the applications of gravity data is geoid modelling. These data are acquired using either satellite observation or gravity 
meters known as gravimeters. The gravity data, whether acquired by satellite means or gravimeter, they are normally reduced to 
the geoid. The geoid according to Eteje et al (2018), is the surface which coincides with that surface to which the oceans would 
conform over the entire earth, if free to adjust to the combined effects of the earth's mass attraction (gravitation) and the 
centrifugal force of the Earth's rotation. It is a surface adopted as a reference for vertical coordinate system. Reducing the gravity 
data/observation to the geoid implies that the observation is astronomical. To apply these observations for geoid model 
determination, their corresponding theoretical gravity of the observed points on a specified ellipsoid is required. The theoretical 
gravity is a latitude dependent component on the ellipsoid. Thus, it is geodetically as well as geometrically obtained. The 
difference between the observed/absolute gravity and the theoretical gravity at a point is known as the gravity anomaly. 

The gravity anomalies of observed points can be either global or local depending on the ellipsoid on which they were obtained. If 
the theoretical gravity are computed on the WGS 84 (global) ellipsoid, the theoretical gravity are said to be global. On the other 
hand, if they are computed on the local ellipsoid adopted for geodetic computation in a country or region, the theoretical gravity 
are said to be local. The computation of the global gravity anomalies requires the application of the International Gravity 
Formula. The International Gravity Formula as given by Aziz et al (2010) is 


St 


f 

9.7803267714 
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1 + 0.00193185138639 sin 2 <p ' 
(1 -0.00669437999013 sin 2 ^) 1/2 , 
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( 1 ) 


There have been no derived/specified models for computation of normal or theoretical gravity on various local ellipsoids in the 
world. Precise practical local geoid determination requires the theoretical gravity to be computed on the specified ellipsoid 
adopted for geodetic computation in the area/region of study. Computing the normal gravity on the local ellipsoid enables the 
gravity anomalies of observed points to be locally determined. Clarke 1880 is the Nigeria Minna Datum ellipsoid. Several local 
gravimetric geoid models have been determined in various parts of the country, but the gravity anomalies of selected/observed 
points were computed using the International Gravity Formula which is global. These were done as a result of the absence of a 
local theoretical gravity model. Consequently, this paper derives a local theoretical gravity model on the Clarke 1880 ellipsoid for 
practical local geoid model determination. 
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The derivation of the local theoretical gravity model on the Clarke 1880 ellipsoid involves the computation of some parameters 
such as semi-minor axis, b , first eccentricity squared, e , second eccentricity, e , acceleration rate, m, q o , q o , theoretical 

gravity at the equator, g e and theoretical gravity at the pole, g p using the Clarke 1880 ellipsoid parameter such as such as semi¬ 
major axis and flattening, and some other constants such as angular velocity and gravitational constant. 

Clarke 1880 Ellipsoid 

The Clarke 1880 ellipsoid is a local ellipsoid adopted for geodetic computation in Nigeria. The semi-major axis a and flattening,/ 
of the Clarke 1880 ellipsoid as given by Eteje et al (2018) are respectively 6378249.145m and 1/293.465. 

The Gravity Formula 

The general gravity formula known as the Somiglinana's closed formula for normal gravity was used for the derivation of the 
International Gravity Formula on the WGS 84 ellipsoid. Moritz (1980) gave the Somiglinana's closed formula for normal gravity 
as 


7 = 8t 


ag e cos 2 <p + bg p sin 2 <p 
yl a 2 cos 2 ~<p + ~b 2 sin 2 ~(p 


( 2 ) 


Where, y — g T — point theoretical gravity 

g e - theoretical gravity at the equator 
g p - theoretical gravity at the pole 
(p = observed point latitude 
a = semi-major axis of specified ellipsoid 
b = semi-minor axis of specified ellipsoid 

The more convenient form of equation (2) for numerical computation is given by Moritz (1980) as 


= (1 + ^sin 2 cp) 

8t ~ {\-e 2 sin 2 (p) lil 

with k given by Madden (2006) as 

bg 
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k = 


-1 
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ag e 


Where, e is the eccentricity squared of the specified ellipsoid. 

The computation of k requires the computation of the theoretical gravity at the equator and the pole. The theoretical gravity at the 
equator and the pole are respectively computed as (Moritz, 1980) 
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GM = gravitational constant = 3986005 xl0 8 m 3 s~ 2 
co 2 = angular velocity = 7292115 x 1CT 11 rad s’ 1 

Also, the computation of the acceleration rate, m, q o and q o , necessitates the computation of the semi-minor axis, first 

eccentricity squared and second eccentricity squared of the specified ellipsoid. The semi-minor axis of an ellipsoid is computed as 
(Clynch, 2002) 


/ 


a-b 


=>b = a(l-f ) 

while the first eccentricity squared is given as (Li and Gotze, 2001) 
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The second eccentricity squared according to Moritz (1980) is 


12 


2 7 2 

a - b 
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Where, 


f = flattening of the ellipsoid 
e = first eccentricity 


e' 2 = second eccentricity squared 
The conventional abbreviated series expansion for the gravity formula is (Moritz, 1980): 


r = g T =ge 
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Gravity Anomaly 

The gravity anomaly is the difference between the points gravity reduced to the geoid and the latitude as well as the normal 
gravity computed on a specified ellipsoid and corrected for free air and effect of rock. The points gravity anomalies are grouped 

into two: the free air gravity anomalies and the Bouguer gravity anomalies. The free air gravity anomaly, A g fa is computed as 

(Aziz et al, 2010): 


^8 fa 8 Obs g fa 8o (17) 

Where, A g fa is free-air anomaly, g obs is observed gravity and g o is theoretical gravity and g fa = 0.3086[mGa l/m]H. 
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The Bouguer gravity anomaly, A g B is given by Vermeer (2016) as 


Ag B = Ag fa - 0.11 19H 


(18) 


The orthometric heights of the points can be obtained from the Digital Terrain Model of the area. 

Gravimetric Geoid 

The gravimetric method of geoid determination can be carried out by the well-known Stokes-integral, equation (19) and the use of 
accurately determined absolute gravity data (Heiskanen and Moritz, 1967, Eteje, 2015 and Eteje et al, 2018). 


N = 


A-JJassWo- 


(19) 


Where, N is geoid undulation, A g is gravity anomaly, S (l// ) is stokes function, y is normal gravity on the reference ellipsoid 
and R is mean radius of the earth. 


DERIVATION OF THE CLARKE 1880 THEORETICAL GRAVITY MODEL 

The derivation of the Clarke 1880 theoretical gravity model involves the computation of the following: semi-minor axis, b, first 
eccentricity squared, e 2 , second eccentricity squared, e' 2 , acceleration rate, m, q o , q o , theoretical gravity at the equator, g e , 

theoretical gravity at the pole, g p , k , and the presentation of the model in its final forms. Figure 1 shows the flow chart of the 
computation procedures. 



Fig. 1: Flow Chart of the Computation Procedures 


Using the given semi-major axis and flattening of the Clarke 1880 ellipsoid, the angular velocity and gravitational constant, the 
computations were carried out as follow: 

Computation of Semi-minor Axis, b 


www.scm.org 

© 2019, Scientific Research Journal 
http://dx.doi.org/10.31364/SCIRJ/v7.i2.2019.P0219612 
























































Scientific Research Journal (SCIRJ), Volume VII, Issue 2, February 2019 
ISSN 2201-2796 


16 


The semi-minor axis was computed using equation (11) as 




b = 6378249.145 


1 


1 




293.465 


6356514.869m 


Computation of First Eccentricity Squared, e 2 

Using equation (12), the first eccentricity squared was computed as 

, 6378249.145 2 - 6356514.869 2 


6378249.145 

Computation of Second Eccentricity Squared, e' 2 

The second eccentricity squared was computed with equation (13) as 


0.00680351145465242m 


e 2 = 


6378249.145 2 - 6356514.869 2 


= 0.006850115m 


6356514.869 

Computation of Acceleration Rate, m 

The acceleration rate, m was computed using equation (9) as 

(7292115 x 10' 11 ) 2 x6378249.145 2 x 6356514.869 „ 

m = ------= 0.003449778 


3986005xl0 5 


Computation of q o 

With equation (7), q o was computed as 
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q'o = 3x 
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v 0.006850115 J 
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0.082765427 


■ x arctan 0.082765427 


-1 


Computation of q o 


Using equation (8), q o was computed as 


2 q, 


q. 


= 1 +- 

t 0.006850115 J 

0.000150304 

2 

= 0.000075152 


x arctan 0.082765427 - 


3 

0.082765427 


0.000150304 


It is to be noted here that the value of arctan 0.082765427 was converted to radian before it was used. 

Computation of Theoretical Gravity at the Equator, g e 

The theoretical gravity at the equator on the Clarke 1880 ellipsoid, g e was computed using equation (5) as 
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8e = 


3986005xl0 8 


6378249.145x6356514.869 


1-0.003449778 


A 


0.003449778 0.082765427 x 0.002732032 

6 X 0.000075152 j 


= 9.780519381ms 


-2 


Computation of Theoretical Gravity at the Pole, g ; 


With equation (6), the theoretical gravity at the pole, g was computed as 


3986005 xlO 8 


6378249.145" 


1 + 


0.003449778 0.082765427 x 0.002732032 


V 


0.000075152 


= 9.831842301ms" 


Computation of k 


By substituting the obtained values of g e and g p into equation (4), k was computed as 


k - 


6356514.869x9.831842301 

6378249.145x9.780519381 


-1= 0.00182202113732433 


Final Forms of the Clarke 1880 Ellipsoid Theoretical Gravity Model, Y aarke i SS0 = g 


T Clarke 1880 


2 

Also, by substituting the values of g e , k, and e into equation (3), the final form of the Clarke 1880 ellipsoid theoretical gravity 
model A was obtained as 


y Clarke 1 


1880(A) St Clarke 1880 (A) 


= 9.780519381 


(l + 0.00182202113732433sin 2 </>) 
(1 -0.00680351145465242 sin 2 tp) 1 ' 2 


ms 


-2 


Again, using the series expansion formula given in equations (14) to (16), the series expansion form of the model (model B) was 
obtained as follow: 

/* was computed using equation (15) as 

= 9.831842301-9.780M9381 

9.780519381 

while f 4 was computed with equation (16) as 

f, = -x 0.00340756137869933 2 +-0.00340756137869933x0.003449778 
2 2 

= 0.00003519406296951 


Substituting the computed values of /* and / 4 into equation (14), we obtained 

/ = g T = 9.780519381^1 + 0.00524746365716546 sin 2 <p - f x 0.00003519406296951 sin 2 2 <p 
= 9.780519381 (l + 0.00524746365716546sin 2 -0.0000087985157424sin 2 2^)ms’ 2 

Therefore, the series expansion form (model B) of the Clarke 1880 ellipsoid theoretical gravity model is 


ms 
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yciarke 1&&0 (b) = &t clarke 1880 (b) = 9.78051938 l(l + 0.0052474636sin 2 cp - 0.00000B7985sin 2 2^)ms- 2 

Comparison of the Models 

Table 1 and figures 2 and 3 show the comparison of the two model forms. The two model forms, that is, Model A obtained with 
equation (3) and model B obtained with equation (14) as well as using the series expansion formula were compared to determine 
the suitability and agreement of the model forms. The theoretical gravity of six points were computed with their latitudes using 
model A and model B. The computed normal gravity of the points using the two model forms were compared to obtain their 
differences. It can be seen from table 1 that the normal gravity of the six points obtained from model A are identical with those 
obtained from modal B which implies that there is no difference between the computation results of the two model forms. Thus, 
any of the model forms can be applied where Clarke 1880 ellipsoid is adopted as a reference surface for geodetic computation. It 
can also be seen from figures 2 and 3 that the plotted normal gravity from the two model forms have the same shape as well as 
their differences which also implies the agreement and suitability of the two model forms. 


Table 1: Comparison of the Models 


Latitude 

Model A (Equation (3)) 

Model B (Equation (14)) 

Difference 

6.331514607 

9.781140782 

9.781139430 

0.000001352 

6.301664360 

9.781134960 

9.781133620 

0.000001339 

6.287544686 

9.781132215 

9.781130881 

0.000001334 

6.257952119 

9.781126483 

9.781125161 

0.000001321 

6.237847990 

9.781122603 

9.781121290 

0.000001313 

6.214210013 

9.781118058 

9.781116755 

0.000001303 



Model A (Equation (3)) 
Model B (Equation (14)) 


Fig. 2: Plot of the Computed Normal gravity 
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CONCLUSION 

Seeing that the gravimetric method of geoid model determination requires the use of gravity anomalies of selected points, the 
normal gravity of points of known absolute or observed gravity are required to be computed. For precise practical local geoid 
model determination, the normal gravity of the selected points are computed on the local ellipsoid adopted for geodetic 
computation in the region/area of study. Since Clarke 1880 is a local ellipsoid adopted for geodetic computation in some countries 
such as Nigeria, this paper has derived a local theoretical gravity model on the Clarke 1880 ellipsoid for practical local geoid 
model determination. The derived model was presented in two forms: model A and model B. The two model forms were 
compared to determine their suitability as well as agreement. The comparison results show that the two forms of the derived 
model agree. Thus, any of the model forms can be applied for accurate normal gravity computation in regions/countries where 
Clarke 1880 ellipsoid is adopted as a reference surface for geodetic computation. 
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